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I.  INTRODUCTION 


The  purpose  of  this  report  Is  to  develop  some  general  numerical  methods 
for  the  problem  of  antenna  pattern  synthesis,  with  and  without  constraints. 

The  principles  of  pattern  synthesis  are  well  known,  and  many  specific  problems 
have  been  considered  [1].  A  good  discussion  of  the  general  theory  and  of  the 
difficulties  encountered  has  been  given  by  Deschamps  and  Cabayan  [2].  To  over¬ 
come  problems  concerning  the  stability  and  sensitivity  of  the  solution,  they 
propose  to  use  regularization  methods  [3].  In  Hilbert  space,  this  method  is 
equivalent  to  obtaining  a  least-iquares  solution  with  a  constraint  on  the 
source  norm  [4,5].  It  is  closely  related  to  obtaining  a  least-squares  solution 
with  a  constraint  on  the  quality  factor,  as  considered  in  Section  VI  of  this 
report. 

Most  of  the  methods  developed  for  pattern  synthesis  assume  that  the 
radiation  field  is  specified  in  both  magnitude  and  phase.  In  many  cases  only 
the  magnitude  is  of  Interest,  and  the  phase  is  left  unspecified.  This  is  a 
special  case  of  the  so-called  mixed  problems  of  antenna  synthesis,  considered 
by  Bakhrakh  and  Troytskiy  [6].  An  application  of  field  magnitude  pattern 

[1]  For  example,  see  R.  E.  Collin  and  F.  J.  Zucker,  "Antenna  Theory,"  Part  1, 
McGraw-Hill  Book  Co.,  New  York,  1969,  Chap.  7. 

[2]  G.  A.  Deschamps  and  H.  S.  Cabayan,  "Antenna  Synthesis  and  Solution  of 
Inverse  Problems  by  Regularization  Methods,"  IEEE  Trans. .  vol.  AP-20, 

No.  3,  May  1972,  pp.  268-274. 

[3]  A.  N.  Tihonov,  "Solution  of  Incorrectly  Formulated  Problems  and  the 
Regularization  Method,"  Soviet  Mathematics,  vol.  4,  July-December 
1963,  pp.  1035-1038. 

[4]  V.  I.  Popovkin  and  V.  I.  Yelumeyev,  "Optimization  and  Systematization 
of  Solutions  to  Antenna  Synthesis  Problems,"  Radio  Engineering  and 
Electronic  Physics,  vol.  13,  No.  5,  1968,  pp.  682-686. 

[5]  V.  I.  Popovkin,  G.  I.  Shcherbakov,  V.  I.  Yelumeyev,  "Optimum  Solutions 

of  Problems  in  Antenna  Synthesis  Theory,"  Radio  Engineering  and  Electronic 
Physics,  vol.  14,  No.  7,  1969,  pp.  1025-1030. 

[6]  L.  D.  Bakhrakh  and  V.  I.  Troytskiy,  "Mixed  Problems  of  Antenna  Synthesis," 
Radio  Engineering  and  Electronic  Physics,  vol.  12,  No.  3,  March  1967, 

pp.  404-414. 
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synthesis  to  a  specific  problem  has  been  published  by  Choni  [7].  In  Section  IV 
we  give  a  discussion  of  field  magnitude  pattern  synthesis  in  terms  of  numerical 
methods. 

The  procedures  developed  in  this  report  are  general,  being  applicable  to 
any  antenna  system  which  can  be  accurately  analyzed.  However,  to  keep  the 
examples  and  computer  programs  simple,  the  theory  is  applied  explicitly  only 
to  point  sources  arbitrarily  distributed  in  a  plane,  and  to  pattern  synthesis 
in  this  plane.  The  extension  to  arbitrary  N-port  systems  in  three-dimensional 
space,  and  to  pattern  synthesis  over  the  entire  radiation  sphere,  is  straight¬ 
forward  but  tedious  in  detail.  A  special  case  of  three-dimensional  field  magni¬ 
tude  pattern  synthesis  applied  to  wire  scatterers  is  given  in  reference  [8]. 

II.  ANTENNA  SYSTEMS  AND  VECTOR  SPACES 


The  general  relationship  between  the  source  f  of  a  radiating  system  and 
the  field  g  it  produces  on  the  radiation  sphere  can  be  symbolized  by 

Tf  -  g  (1) 

where  T  is  a  linear  operator.  In  general,  both  f  and  g  may  be  infinite 
dimensional  and  T  1  may  not  be  bounded.  The  pattern  synthesis  problem  can  be 
summarized  as  follows:  Given  a  desired  field  gQ  (specified  completely  or 
partially),  we  wish  to  determine  a  source  f  (constrained  or  unconstrained) 
whose  field  g  approximates  gQ  in  some  acceptable  manner.  The  approach  taken 
by  Deschamps  and  Cabayan  [2]  is  to  treat  the  problem  in  general  in  function 
spaces,  and  then  numerically  evaluate  the  final  formulas  for  specific  appli¬ 
cations.  Our  approach  is  to  discretize  the  problem  at  the  beginning,  and  then 
treat  the  problem  in  finite-dimensional  vector  spaces. 


[7]  Y.  I.  Choni,  "Synthesis  of  an  Antenna  According  to  a  Given  Amplitude 
Radiation  Pattern,"  Radio  Engineering  and  Electronic  Physics,  vol.  16, 

No.  5,  May  1971,  pp.  770-778. 

[8]  R.  F.  Harrington  and  J.  R.  Mautz,  "Synthesis  of  Loaded  N-port  Scatterers," 
AFCRL-72-0665,  Scientific  Report  No.  17  on  Contract  No.  F19628-68-C-0180 
between  Syracuse  University  and  Air  Force  Cambridge  Research  Laboratories, 
October  1972. 
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To  render  Che  source  discrete,  we  assume  Chat 


f 


N 

l  f 


n"l 


e 

n  n 


(2) 


where  f  are  constants  and  e  are  basis  elements.  If  the  source  is 
n  n 

continuous,  the  e  are  functions  and  (2)  is  usually  an  approximation  to' 
n 

the  true  source.  If  the  source  is  discrete,  the  e  are  finite-dimensional 

n  * 

vectors,  and  (2)  is  an  exact  relationship.  We  define  the  vector  f  to  be  the 
vector  of  the  components  f^,  that  is 


1 


*f  nVl 


(3) 


For  example,  if  the  source  is  an  array  of  dipoles,  the  f  may  be  the  input 
port  currents,  or  the  input  port  voltages,  or  their  components  in  some 
arbitrary  basis.  We  substitute  (2)  into  (1)  and  evaluate  the  equation  at 

M  points  (&!„»♦„,)»  n>“l,2 . M  on  the  radiation  sphere.  The  result  can  be 

written  as  the  matrix  equation 


where  g  is  the  vector 


[Til  ■  t 


8  -  [gm]Mxl 


(4) 

(5) 


and  [T]  is  the  matrix 


[T] 


^Ten^m^MxN 


(6) 


Here  g  denotes  the  value  of  g  at  the  point  (0  ,  $  ),  and  (Te  )  denotes 
in  m  ro  n  m 

the  pattern  of  f^  evaluated  at  the  point  these  quantities  are 

spatial  vectors,  then  the  above  procedure  applies  to  the  6  and  <j>  components 
of  the  vector  field.  More  generally,  one  could  define  a  set  of  testing 
functions  and  reduce  (1)  to  a  matrix  equation  by  the  method  of  moments  [9]. 


[9]  R.  F.  Harrington,  "Field  Computation  by  Moment  Methods,"  Macmillan  Co., 
New  York,  1968. 
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The  synthesis  problem  cen  now  be  represented  by 

(T]*«g0  (7) 

where  gQ  Is  the  specified  pattern  vector.  Equation  (7)  is  a  matrix  equa¬ 
tion,  and  the  well-known  methods  of  matrix  algebra  can  now  be  brought  to 
bear  on  the  problem.  If  [T]  is  square  and  nonsingular,  the  solution  is 
given  by  the  inverse  to  (7) .  This  solution  will  give  the  correct  value  of 
field  at  the  M  points  (6  )  on  the  radiation  sphere,  but  no  control  of 

tQ  IB 

the  field  is  obtained  between  these  points.  If  [T]  is  of  rank  less  than  N, 
the  problem  is  underdetermined  and  more  than  one  solution  to  (7)  exists.  In 
this  case  additional  constraints  can  be  applied  to  either  the  source  or  the 
field  to  obtain  the  most  desirable  solution.  If  [T]  is  of  rank  greater  than 
N,  then  usually  no  exact  solution  to  (7)  exists,  but  there  will  be  a  uniaue 
least-squares  solution.  Again  additional  constraints  can  be  applied  to  ob¬ 
tain  more  desirable  solutions. 

Let  us  next  define  the  Euclidean  spaces  involved  more  precisely.  For 
source  quantities,  the  inner  product  is  defined  as 

<fi’  V  "  *i[V1?j  (8) 

where  the  tilde  denotes  transpose,  the  asterisk  denotes  complex  conjugate, 
and  [V]  is  a  weight  matrix.  [V]  must  be  positive  definite,  and  for  this 
report  we  consider  it  to  be  the  diagonal  matrix 

'v>  '  !“«  <’> 

with  all  v^  >  0.  In  the  source  space  the  norm  induced  by  the  inner  pro¬ 
duct  is 

Util  .  .  (  J  »  If  |2>1'2  (10) 

n-1 

where  fR  are  the  components  of  1.  The  metric  induced  by  the  norm  is 
d(?^,lj)  ■  ||  f^  -  ||  ■  For  field  quantities  we  define  the  inner  pro¬ 

duct  as 
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(11) 


<81«8j>  -  ^[W]gj 

where  [W]  Is  a  weight  matrix.  It  must  be  positive  definite  and  for  this 
report  we  take  it  to  be  the  diagonal  matrix 

[W]  -  [diag  wm]MxM  (12) 

with  all  wffl  >  0.  In  the  field  space  the  iorm  induced  by  the  inner  pro¬ 
duct  is 

llgll  *  <g*«8>1^2  “  (  l  wjsj2)1^2  (13) 

m-1 

where  g  ere  the  components  of  g.  The  metric  Induced  by  the  norm  is 

d^.gj)  -Ilg^jH- 

To  compare  various  pattern  synthesis  procedures,  we  define  two 
figures  of  merit.  The  first  is  the  normalized  synthesis  error 


E 


M 

l  w 

m-1 


g  -g  I' 
m'°m  om1 


M 


l  wm  l» 


m-1 


(14) 


where  g  is  the  synthesized  pattern  and  gQ  is  the  desired  pattern.  The 
second  is  the  quality  factor  [10] 


Q 


1a|2 

"  "A'2 

m-1 


(15) 


The  multiplier  M  is  introduced  into  (15)  to  make  the  Q  relatively  in¬ 
sensitive  to  the  number  of  field  points  chosen. 


[10]  Y.  T.  Lo,  S.  W.  Lee,  and  Q.  H.  Lee,  "Optimization  of  Directivity  and 
Signal-to-Noise  Ratio  of  an  Arbitrary  Antenna  Array,  "  Proc.  IEEE, 
vol.  54,  August  1966,  pp.  1033-1045. 
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III.  PATTERN  SYNTHESIS  WITH  PHASE  SPECIFIED 


The  most  commonly  used  method  of  pattern  synthesis  involves  specifying 
the  field  pattern  in  both  magnitude  and  phase.  For  our  method,  this  involves 
specifying  the  magnitude  and  phase  of  gQ  at  M  points  on  the  radiation  sphere. 
Hence,  the  starting  point  is  (7)  with  gQ  known.  We  normally  take  more  inde¬ 
pendent  equations  than  unknowns,  in  which  case  the  least-squares  solution  to 
(7)  is 


1  -  [T*WT]_1[Y*W]go  (16) 

This  is  obtained  in  the  usual  way  by  minimizing  the  quantity  ||  [Tjl-g^H2. 
The  derivation  of  (16)  follows  (31)  to  (34)  of  Section  IV  with  o-O. 
If  unweighted  inner  products  are  used,  then  [W]  is  the  identity  matrix. 


For  examples,  the  equations  are  specialized  to  point  sources  arbitrarily 

distributed  in  a  plane  and  to  radiation  patterns  in  this  same  plane.  Figure 

1  illustrates  the  general  problem.  The  point  sources  are  located  at  positions 

(x  ,y  )  and  have  excitations  f  .  The  radiation  pattern  at  the  angle  <j>  is  then 
n  n  n 

given  by 


N  Jk(x  cos  ♦  +  y  sin  if) 

Tf-  J  f«  n  n  (17) 

Li  n 
n-1 

where  k  »  2u/A  is  the  wavenumber.  Choosing  a  desired  pattern  gQ(<f),  and 

setting  Tf  *  go  at  H  points  we  obtain  [T]l  *sgQ.  Here  ?  is  the  column 

vector  of  the  f  ,  g  is  the  column  vector  of  the  g  (f  ) ,  and  [T]  ie  the 
no  o  m 

matrix  with  elements 


T 

mn 


ik(x  cos  $  +  y  sin  $  ) 

n  m  n  m 


(18) 


For  simplicity  we  take  the  unweighted  inner  product,  in  which  case  [W]  is 
the  identity  matrix.  The  solution  for  the  sources  f  is  then  given  by  (16), 
and  the  synthesized  pattern  by  (17).  This  solution  is  valid  so  long  as  M 
is  sufficiently  large  to  give  at  least  N  independent  equations  in  the  N  un¬ 
knowns  t  . 

n 
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For  some  numerical  results,  consider  10  sources  approximately 

equispaced  on  one-half  of  a  2:1  ellipse,  as  shown  in  Fig.  2a.  For  the 

radiation  pattern,  take  a  rectangle  in  the  first  quadrant  of  the  polar 

pattern.  (This  is  the  so-called  cosecant  $  pattern.)  36  field  points 

$  are  chosen  every  10s  in  $,  starting  at  $  »  5°.  These  are  shown  by 
m 

crosses  on  Fig.  2b.  For  the  synthesis  we  take  four  cases:  (a)  origirj 
at  the  center  of  the  ellipse  and  the  field  real,  (b)  origin  one-half 
the  distance  from  the  center  to  the  end  point  of  the  semimajor  axis  and 
the  field  real,  (c)  origin  at  the  end  point  of  the  semimajor  axis  and  the 
field  real,  and  (d)  origin  at  the  center  of  the  ellipse  and  the  field  phase 
alternating  from  0  to  180s  between  adjacent  field  points.  (This  last  case 
represents  the  worst  choice  that  can  be  made.)  Figure  3  shows  the  results 
when  the  sources  are  separated  d  «  A/4,  Fig.  4  when  d  -  A/2,  and  Fig.  5  when 
d  -  A.  The  corresponding  normalized  synthesis  errors  E  and  quality  factors  0 
associated  with  each  synthesis  result  are  given  below  each  pattern.  The 
source  excitations  are  listed  by  magnitude  and  phase  (in  degrees)  in  Tables  1 
to  3,  normalized  so  that  the  maximum  excitation  is  unity . 

Some  general  observations  are  as  follows:  When  the  field  is  chosen  to 
have  alternating  phase  between  adjacent  points  (cases  d),  the  pattern  syn¬ 
thesis  is  poor  in  each  case.  When  the  field  is  chosen  to  be  in  phase  at  all 
points,  with  respect  to  a  coordinate  origin  in  the  vicinity  of  the  sources, 
the  synthesis  is  good  when  the  sources  are  separated  by  d  “  A/4  and  A/2,  and 
poor  when  separated  by  d  *  A.  The  Q  is  low  when  the  sources  are  separated  by 
d  -  A/2  and  A,  but  relatively  high  when  d  -  A/4.  These  properties  are  in  agree¬ 
ment  with  what  we  would  expect  based  on  past  experience  with  antenna  synthesis 
problems. 
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(c)  E  -  0.223,  Q  =  11.1  (d)  F.  -  0.957,  Q  *  20.5 


Fig.  3.  Antenna  pattern  synthesis,  phase  specified,  d  -  1/4. 
Specified  pattern  shown  by  crosses,  synthesized 
pattern  shown  solid. 
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Fig.  A.  Antenna  pattern  synthesis,  phase  specified,  d  -  \/2. 
Specified  pattern  shown  by  crosses,  synthesized 
pattern  shown  solid. 
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Fig.  S.  Antenna  pattern  synthesis,  phase  specified,  d  *  X. 
Specified  pattern  shown  by  crosses,  synthesized 
pattern  shown  solid. 
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Table  1.  Element  excitations  for  the  synthesized  patterns  of  Fig.  3. 


Element 

Figure  3a 
mag.  phase 

Figure  3b 
mag.  phase 

Figure  3c 
mag.  phase 

Figure  3d 
mag.  phase 

fl 

0.312 

-38.3* 

0.164 

84.2s 

0.254 

-29.2s 

0.441 

-26.7s 

f2 

0.762 

173.9s 

0.627 

-73.0s 

0.657 

159.2s 

0.858 

167.0s 

f3 

1.000 

0 

0.810 

126,4s 

1.000 

0 

1.000 

0 

f4 

0.903 

-147.4s 

0.761 

-27.8s 

0.955 

-147.1s 

0.680 

178.1s 

f5 

0.502 

32.6s 

0.482 

155.4s 

0.462 

55.8s 

0.394 

-10.3s 

f6 

0.605 

-122.3s 

0.707 

-5.4s 

0.529 

-137.2s 

0.203 

103.0s 

f7 

0.507 

37.1s 

0.663 

178.9s 

0.580 

43.1s 

0.363 

-138.6s 

f8 

0.717 

-126.8s 

1.000 

0 

0.868 

-147.8° 

0.669 

25.6s 

f9 

0.647 

16.4s 

0.664 

164.1s 

0.678 

17.3s 

0.612 

179.8s 

f10 

0.300 

149.8s 

0.300 

-34.3s 

0.331 

164.9s 

0.292 

-17.6s 

Table  2.  Element  excltatlone  for  the  ayntheaized  pattema  of  Fig.  4 


Element 
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Figure  4a 
mag.  phase 

Figure  4b 
mag.  phase 

Figure  4c 
mag.  phase 

Figure  4d 
mag.  phase 

0.341 

60.9° 

0.525 

-115.4s 

0.383 

-90.6' 

0.713 

-133.0° 

0.315 

-66.7° 

0.162 

-11.6s 

0.340 

77.5s 

0.502 

46.8s 

0.378 

22.5s 

0.231 

64.0° 

0.297 

-18.6s 

0.207 

167.3s 

0.488 

48.8s 

0.248 

132.9s 

0.392 

167.4s 

0.623 

65.6s 

0.767 

161.2s 

0.66' 

-45.6s 

1.000 

0 

1.000 

0 

0.989 

-179.0s 

0.939 

-42.1s 

0.847 

-37.3s 

0.783 

8.6s 

0.960 

-112.7° 

1.000 

0 

0.427 

36.8s 

0.510 

57.1s 

1.000 

0 

0.427 

66.5s 

0.365 

136.7s 

0.039 

-74.5s 

0.630 

94.6s 

0.305 

-25.9s 

0.341 

13.8s 

0.547 

35.2s 

0.057 

-75.7s 

0.102 

130.5s 

0.240 

175.7s 

0.734 

-123.9s 
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Table  3.  Element  excitations  for  the  synthesized  patterns  of  Fig.  5. 


Element 

Figure  3a 
mag.  phase 

Figure  5b 
mag .  phe  se 

Figure  5c 
mag.  phase 

Figure  5d 
mag.  phase 

fl 

— 

0.062 

-173.5* 

0.252 

50.6° 

0.448 

-134.0° 

0.619 

17.3* 

f2 

0.063 

18.6* 

0.384 

-129.6* 

0.284 

-108.1° 

0.663 

21.1’ 

f3 

0.510 

102.9* 

0.467 

-4.4° 

0.508 

-125.1° 

0.134 

10.1' 

f4 

0.377 

130.3° 

0.458 

-22.2° 

0.379 

-82.6“ 

]  .000 

0 

f5 

0.497 

-94.5* 

0.551 

120.4° 

1.000 

0 

0.284 

1.2’ 

f6 

0.939 

153.6° 

1.000 

0 

0.815 

-150.1° 

0.700 

86.1° 

f  7 

1.000 

0 

0.802 

168.8° 

0.366 

58.3° 

0.150 

-98.5° 

f8 

0.848 

-107.3° 

0.442 

148.5* 

0.481 

32.7° 

0.247 

26.1° 

f9 

0.285 

34.9° 

0.129 

28.5° 

0.151 

-64.4° 

0.781 

4.3’ 

f10 

0.120 

-168.5° 

0.068 

118.9° 

0.126 

147.4° 

0.499 

? .  3° 

IV.  FIELD  MAGNITUDE  PATTERN  SYNTHESIS 


We  next  consider  the  problem  of  synthesizing  a  pattern  in  magnitude 
only.  Let  h  -  |gQ|  be  a  desired  field  magnitude,  and  form  the  vector  h  by 
specifying  its  value  h  at  M  points  on  the  radiation  sphere.  Again  we  con- 
aider  the  source  to  be  discretized  and  represented  by  the  vector  f.  It  is 
desired  to  find  the  source  t  for  which  the  pattern  error 


e  -  ||  |[T]?|  -  £||2  (13) 

is  minimum.  In  terms  of  components,  (19)  becomes 


where  the  w  are  weight  factors.  To  circumvent  the  troublesome  inner 
m 

magnitude  operation  in  (20),  we  first  consider  the  more  general  function 


e(f,1) 


(21) 


This  is  the  error  function  used  when  the  pattern  is  specified  in  both  magnitude 

h  and  in  phase  0  .  Hence,  for  6  fixed,  the  f  for  minimum  e  are  given  by 
in  id  m  n 

(16).  For  f^  fixed,  the  minimum  e  is  obtained  when  both  terms  within  the  mag¬ 
nitude  signs  of  (21)  are  in  phase,  that  is,  when 


(22) 


Because  (21)  is  more  general  that  (20),  its  minimum  is  less  than  or  equal  to 
that  of  (20).  But  under  condition  (22),  the  e  of  (21)  is  equal  to  that  of 
(20).  Therefore  (20)  and  (21)  have  the  same  minimum. 
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An  Iterative  procedure  for  minimizing  (21)  proceeds  as  follows: 


1.  Assume  starting  values  for  8^,  S^t**** 

2.  Keep  the  B  fixed  and  calculate  the  f  which  minimize  e  using  (16). 

in  n 

3.  Keep  the  f  fixed  and  calculate  the  B  which  minimize  e  using  (22). 

n  in 

4.  Go  to  step  2. 


This  procedure  eventually  converges  because  steps  2  and  3  cannot  increase  e. 
While  the  procedure  obtains  absolute  minima  in  the  ?  space  and  in  the  8  space, 
it  does  not  necessarily  obtain  the  absolute  minimum  in  the  catenated  space 
a.h.  Hence,  the  procedure  converges  to  a  stationary  point,  usually  a  local 
minimum,  which  may  or  may  not  be  the  global  minimum.  An  alternative  procedure 
for  minimizing  (21)  is  given  in  the  Appendix. 


For  numerical  results,  we  consider  the  same  example  as  used  in  the  pre¬ 
ceding  section.  Hence,  the  array  is  illustrated  by  Fig.  2a,  and  the  radiation 
pattern  by  Fig.  2b.  The  same  three  cases  of  element  separation,  d  -  A/4,  X/2, 
and  A,  are  used.  Two  starting  points  were  chosen  for  the  iterative  procedure: 
(a)  origin  at  the  end  point  of  the  semimajor  axis  and  the  field  real,  and  (b) 
origin  at  the  center  of  the  ellipse  and  the  field  phase  alternating  between  0 
and  180s  between  adjacent  field  points.  These  starting  points  correspond  to 
cases  (c)  and  (d)  of  the  previous  section.  The  final  results  of  the  magnitude 
synthesis  procedure  are  shown  in  Fig.  6  for  d  »  A/4,  Fig.  7  for  d  »  A/2,  and 
Fig.  8  for  d  «  A.  The  normalized  synthesis  errors  E  and  quality  factors  Q  for 
each  result  are  given  below  each  pattern.  The  source  excitations  are  listed 
by  magnitude  and  phase  in  Table  4,  normalized  so  that  the  maximum  excitation 
is  unity. 
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Fig.  6.  Field  magnitude  pattern  synthesis,  d  »  1/4. 

(a)  Starting  field  taken  real,  (b)  Starting 
field  alternating  0  to  180°  in  phase  between 
adjacent  field  points. 
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(b)  E  -  0.425,  Q  =  1.09 


Fig.  7.  £ield  magnitude  pattern  synthesis,  d  =  V/2. 

(a)  Starting  field  taken  real,  (b)  Starting 
field  alternating  0  to  180°  in  phase  between 
adjacent  field  points. 
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(b)  E  -  0.517,  Q  -  0.943 


Fig.  8.  Field  magnitude  pattern  synthesis,  d  -  A. 

(a)  Starting  field  taken  real,  (b)  Starting 
field  alternating  0  to  180°  in  phase  between 
adjacent  field  points. 
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Table  4.  Element  excitation*  for  the  synthesized  patterns  of  Figs.  6,  7,  8. 

Ele-j  Figure  6a  |  Figure  6b  |  Figure  7a  |  Figure  7b  |  Figure  8a  |  Figure  8b 

A 


fl 

0.264  125.0 

f2 

0.626 

-40.5 

f3 

0.985 

153.2 

f4 

1.000 

0.0 

f5 

0.636  -159.1 

f6 

0.553 

24.7 

f7 

0.570  -160.3 

f8 

0.664 

10.3 

f9 

0.527 

172.9 

f10 

0.234  -42.8 

0.285  -177,4 


The  following  are  some  observations  concerning  the  improvement  of 
magnitude  pattern  synthesis  over  ordinary  pattern  synthesis.  For  the  cases 
d  «  A/4  and  d  «  A/2,  when  the  starting  field  is  taken  real  (equiphase) ,  the 
improvement  is  relatively  small.  (Compare  Figs.  3c  to  6a,  and  4c  to  7a.) 
For  the  cases  d  -  A/4  and  d  -  A/2,  when  the  starting  field  alternates  0  to 
180°  between  adjacent  field  points,  the  Improvement  is  large.  (Compare 
Figs.  3d  to  6b  and  4d  to  7b.)  For  the  case  d  »  A/4  the  two  final  patterns 
are  the  same  (Figs.  6a  and  6b),  but  for  the  case  d  »  A/2  the  two  final 
patterns  are  different  (Figs.  7a  and  7b).  For  the  case  d  »  A,  the  improve¬ 
ment  over  the  starting  pattern  is  larger  in  both  cases.  (Compare  Figs.  5c 
to  8a  and  5d  to  8b.)  However,  in  no  case  is  the  final  pattern  very  good 
when  d  **  A.  Apparently  the  source  separation  is  too  great  for  good  pattern 
synthesis.  Finally,  note  that  for  the  case  d  ■  A  the  two  final  patterns  are 
the  same  (Figs.  8a  and  8b). 

For  future  reference,  the  source  norm  squared,  the  normalized  syn¬ 
thesis  error,  and  the  Q  for  each  result  are  tabulated  in  Table  5  for  the 
case  d  ■  A/4.  The  five  rows  correspond  to  the  synthesized  patterns  of 
Figs.  3a,  b,  c,  d,  and  Fig.  6.  When  the  source  separation  is  large  the 
various  quantities  of  Table  5  are  less  sensitive  to  the  type  of  synthesis 
used,  and  hence  we  do  not  tabulate  them  for  d  ■  A/2  and  d  =  A. 

Table  5.  Source  norm  squared,  normalized  synthesis  error,  and  quality 
factor  for  unconstrained  pattern  synthesis,  d  »  A/4. 


Ilf  IF 

E 

Q 

coordinate  origin  a 

13.37 

0.312 

12.6 

coordinate  origin  b 

17.14 

0.307 

16.0 

coordinate  origi  :  c 

13.35 

0.223 

11.1 

field  phase  alternating 

1.38 

C.957 

20.5 

magnitude  synthesis 

27.50 

0.172 

21.5 
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V.  PATTERN  SYNTHESIS  WITH  CONSTRAINED  SOURCE  NORM 


The  source  norm  is  closely  related  to  near-field  quantities, 
such  as  power  losses  in  an  antenna  structure  or  energy  storage.  For  this 
reason  it  is  often  desirable  to  limit  the  source  norm,  especially  when  the 
sources  are  close  together  or  continuously  distributed.  Hence,  we  consider 
the  problem  of  minimizing 

e  -  ||  [T]?  -  gjl2  (23) 

subject  to  the  constraint 

||f  ||2  1  C  (24) 

where  C  is  a  positive  constant  to  be  chosen.  This  constrained  minimi¬ 
zation  can  be  accomplished  by  forming  the  Lagrangian 

J  -  II  [Tl?  -  IJI2  +  a  ||1||2  (25) 

where  a  is  a  Lagrange  multiplier.  If  we  can  obtain  the  source  function  f 
which  minimizes  ,J  with  respect  to  f  and  are  able  to  find  a  >  0  such  that 

II  ? H  2  =  C  (26) 

then  any  other  source  function  which  satisfies  the  constraint  (24)  gives 
at  least  as  large  a  J  and  thus  at  least  as  large  an  e  as  that  provided 
by  ?.  Hence  1  minimizes  e  subjec.  to  the  constraint  (24). 

However,  it  is  not  always  possible  to  find  a  >  0  such  that  the 
f  which  minimizes  J  satisfies  (26).  This  warrants  an  investigation  into 
the  behavior  of  || ? j| 2  and  e  attained  by  the  minimizing  function  f  versus 
a.  If  is  the  minimizing  function  when  a  =  and  ^  is  the  minimizing 
function  when  a  =  a^,  then 
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'l  +  »l  II  ^lH 2  —  c2  +  °1  H?2II2 

(27) 

c2 +  °2  h*2ii!5‘i  +  «2  iiy2 

(28) 

where  or  e 2  is  the  e  of  (23)  with  f  replaced  by  or 

inequalities  (27)  and  (28) ,  we  have 

Adding  the 

(a2  -  «x)(  l|f2ll2-  li^ll  2)  £0 

(29) 

which  shows  that  1 1 ? ||2  is  a  monotone  decreasing  function  of  a.  If  both 
a^  and  a2  are  positive,  (27)  and  (28)  can  be  multiplied  by  a2  and  a^,  respec 
tively,  without  changing  the  sense  of  the  inequalities.  The  resulting 

inequalities  can  now  be  added  to  obtain 

(a2  -  a1)(e2  -£^0 

(30) 

which  shows  that  e  is  a  monotone  increasing  function  of  a  for  positive  a. 
The  smallest  possible  e  is  obtained  when  a  -  0,  corresponding  to  the  uncon¬ 
strained  optimization.  The  smallest  possible  ||?||2  is  obtained  when  a  *■  », 
because  when  u  is  very  large  the  Lagrangian  is  essentially  a  ||1||2  and  thus 
minimization  of  the  Lagrangian  is  the  same  as  minimization  of  ||?||2, 

Now  if  it  is  impossible  to  find  a  >  0  such  that  the  minimizing  ? 

satisfies  (26),  then  C  is  either  larger  than  ||f||2  when  a  »  0,  or  C  is 
less  than  ||?||2  when  a  »  «°.  If  C  is  larger  than  ||1||2  when  a  «  0,  then 

the  desired  ?  is  the  one  which  minimizes  J  when  a  ■  0,  because  this  ? 

gives  the  minimum  possible  e  and  does  not  violate  the  constraint  (24) . 
Actually,  the  constraint  (24)  is  ineffective  whenever  C  is  larger  than 
||  f||2  when  a  ■*  0.  If  C  is  smaller  than  ||?||2  when  a  -  ”,  then  since  ||?||2  at 
a  »  •  is  the  smallest  possible  ||?||2  ,  the  constrained  minimization  problem 
(23),  (24)  has  no  solution  because  it  is  impossible  to  satisfy  the  con¬ 
straint  (24). 
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Consider  the  minimization  of  (25)  for  a  >  0  such  that  (26)  is 
satisfied.  For  a  fixed  a  >  0,  (25)  has  an  absolute  minimum  because 
J  >  0.  The  Lagrangian  J,  being  a  quadratic  form  in  f,  has  one  station¬ 
ary  point.  The  absolute  minimum  of  J  mi-st  occur  at  this  stationary 
point  which  can  be  obtained  by  setting  to  zero  the  variations  of  J  with 
respect  to  both  the  real  and  imaginary  parts  of  f,  or  equivalently  ns 
on  page  192  of  reference  [9,  with  respect  to  ?  and  ?  . 

To  minimize  (25),  expand  it  in  terms  of  inner  products  as 

J  -  UT)f  :  g0)*[W]([T]?  -  1o)  +  af*[V]?  (31) 

where  [V]  and  [W]  are  the  weight  matrices  of  (8)  and  (11),  respectively. 
Further  expanding  (31),  we  have 

aJ*  .  ,•>  /*/k  A/ft  -► 

J  -  f  [T  WT]f  -  go(WT]£  -  f  [T  W]gQ 

+  g*lW]g0  +“?*[V}f  (32) 

r>Sk 

Taking  the  variation  of  this  with  respect  to  f  and  setting  it  equal  to 
zero,  we  obtain 

6J  -  6f* ([T*WT]1  -  [T*W]gQ  +  a[V]f )  -  0  (33) 

•A/ft 

Since  6f  is  arbitrary,  the  vector  in  the  parentheses  must  be  zero,  or 

[T*WT  +  aV]f  »  lT*W]gQ  (34) 

Taking  the  variation  of  (32)  with  respect  to  ?  and  setting  it  equal  to 

zero  we  obtain  the  conjugate  equation  to  (33).  Note  that  a  ■  0  in  (34) _ 

[9]  R.  F.  Harrington,  "Field  Computation  by  Moment  Methods,"  Macmillan  Co., 
New  York,  1968. 
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gives  the  usual  least-squares  solution  (equation  (16)  of  Section 
III). 

Next  we  have  to  determine  a  given  the  constraint  ||?|p  -  C.  For  this, 
we  first  consider  the  eigenvalue  equation 

[T*WT)«1  -  Xi(V]*1  (35) 

Let  the  eigenfunctions  be  normalized  with  weight  [V],  so  the  orthogonality 
relationships  become 


■  6iJ 

(36) 

k  -> 

♦  1LT  WTJ*.,  ■  X1  *ij 

(37) 

a /k 


The  matrix  [T  WT]  is  Hermitian  and  [V]  is  positive  definite,  therefore  the 
Ji^  form  r  complete  set  in  the  f  space  and  we  can  write 


f-  I*A 

i-1  1 


(38) 


where  the  a^  are  constants.  Substituting  (38)  into  (34),  we  have 


N 

I 

i-1 


I  ai[T  WT  +  aV]^  -  [T  W]gQ 


(39) 


Premultiplying  (39)  by  <t>j  >  and  using  the  orthogonality  relationships 
(36)  and  (37),  we  obtain 


.  /v*  *sk  -* 

a1(X1  +  a)  -  4>  1  tT  WjgQ 


(40) 


which  determines  the  a^.  Substituting  these  a A  into  (38),  we  have 


u  ?  -A_; 

i-i  xi  +  a  1 


(41) 
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where  Che  constants  are 


MT  Wk 


(42) 


Next,  substitute  (41)  Into  (23),  again  use  orthogonality  of  the  $  ,  and 
obtain 


N  i 

+  l 

1-1  1 


(43) 


Also,  from  (41)  and  orthogonality,  we  obtain 


ll*IP 


N 


l 

1-1 


(Xi  +  a)7 


(44) 


Now  [V]  is  positive  definite  and  f¥*WT]  Is  at  least  positive  indefinite, 
therefore  all  A^  >_  0.  Hence,  as  expected,  the  error  (43)  is  a  monotone 
increasing  function  of  a  for  a  >_  0  and  source  norm  squared  (44)  is  a  mono¬ 
tone  decreasing  function  of  a  for  a  ^  0. 


Because  of  the  monotone  decreasing  nature  of  (44),  there  is  pre¬ 
cisely  one  a  >  0,  say  a ,  which  satisfies 


»  M2 

c '  Ji  ' F<a> ' 0 


(45) 


This  can  be  computed  using  Newton's  method. 

The  solution  i  which  minimizes  (25)  for  a  >  0  such  that  (26)  is 
satisfied  is  unique  and  is  given  by  (41)  with  a  -  a  .  If  C  is  neither 
too  large  nor  too  small,  this  ?  minimizes  e  subject  to  the  constraint  (24). 
As  mentioned  before,  if  C  is  larger  than  the  norm  squared  of  the  ?  which 
minimizes  e  when  there  is  no  constraint,  then  this  ?  minimizes  e  subject  to 
the  ineffective  constraint  (24).  If  C  is  smaller  than||f|p  at  a  -  ",  then 
it  is  impossible  to  satisfy  the  constraint  (24) . 

For  examples,  the  previously  synthesized  patterns  for  the  case  d  -  A/4 
were  rerun  with  a  constraint  on  thp  source  norm.  The  resulting  synthesized 
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Fig.  9.  Pattern  synthesis  with  constrained  source  norm,  d  ”  X/4. 

Field  phase  is  specified  in  (a),  (b) ,  and  (c).  Field  is 
specified  only  in  magnitude  in  (d). 
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patterns  are  shown  in  Fig.  9.  The  first  three  cases  are  for  the  field  real 
with  coordinate  origin  at  the  points  (a),  (b),  and  (c)  of  Fig.  2a.  In  each 
case  the  constraint  was  ||?||2  ■  4.  The  corresponding  patterns  for  uncon¬ 
strained  synthesis  are  those  of  Figs.  3(a),  (b) ,  and  (c).  Note  that  the  final 
synthesized  patterns  are  not  greatly  different  from  the  unconstrained  results, 
yet  ||?|p  has  been  reduced  from  the  order  of  15  to  4  (see  Table  5).  The  choice 
of  coordinate  origin  (a)  and  the  field  phase  alternating  between  adjacent  field 
points  was  also  run  with  the  constraint  ||?|p  ■  1,  and  the  resulting  pattern 
war  essentially  the  same  as  Fig.  3(d).  Finally,  field  magnitude  pattern  syn¬ 
thesis  with  the  constraint  ||f|p  -  8  was  run  using  each  of  the  above  mentioned 
four  starting  points.  The  final  synthesized  pattern  was  the  same  regardless 
of  the  starting  point,  and  the  result  is  shown  in  Fig.  9(d).  Again  the  syn¬ 
thesized  pattern  is  not  greatly  different  from  the  unconstrained  result,  Fig. 

3,  even  though  ||f||2  has  been  reduced  from  27.5  to  8. 

Table  6  lists  ||f||2  E,  and  Q  for  the  constrained  synthesis  results.  It 
should  be  compared  to  Table  5  for  the  corresponding  unconstrained  results. 

Note  that  the  errors  for  the  constrained  patterns  are  always  as  high  or  higher 
than  those  for  the  unconstrained  patterns,  as  they  must  be.  Note  also  that 
the  Q's  of  the  constrained  patterns  are  always  as  low  or  lower  than  those  for 
the  unconstrained  patterns.  We  show  in  the  next  section  that  this  is  generally 
true . 

Table  6.  Source  norm  squared,  normalized  pattern  error,  and  quality  factor 
for  pattern  synthesis  with  constrained  source  norm. 


csl 

til 

E 

Q 

coordinate  origin  a 

4.00 

0.324 

4.05 

coordinate  origin  b 

4.00 

0.326 

4.09 

coordinate  origin  c 

4.00 

0.235 

3.55 

field  phase  alternating 

1.00 

0.957 

16.15 

magnitude  synthesis 

8.00 

0.191 

6.82 
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VI.  PATTERN  SYNTHESIS  WITH  CONSTRAINED  QUALITY  FACTOR 


It  Is  sometimes  desirable  to  constrain  the  quality  factor  Q  defined 
by  (IS).  According  to  the  general  Lagrange  multiplier  theory  in  section  V, 
the  pattern  synthesis  error  c  can  be  minimized  subject  to  the  constraint 

Q  <  Q0  (46) 

by  minimizing  a  Lagranglan 

J  -  e  +  aQ  (47) 

with  respect  to  the  source  function  ?  and  choosing  a  >  0  such  that 


Q  -  Q_ 


(48) 


at  least  if  Qq  is  neither  too  large  nor  too  small. 


With  (23)  and  (15),  the  Lagranglan  (47)  is  given  in  terms  of  ?  by 

(49) 


„  *  .  M2 


Similar  to  (32),  we  have 


f  [V]f 


J  -  f  [T  WT]f  -  go[WT]f-  f  [T  W]gQ+  g0[W]gQ  +  oM  (50) 

Taking  the  variation  of  (50)  with  respect  to  f  and  setting  it  equal  to 
zero,  we  have 


6J  -  <5f*((l  -  ..  ■  a3...  )[T*WT]1  +  .— ~r~T  (Vjl  -  [T* W]g  )“0  (51) 


[T  ]?|p 


II  [Tjf 


Since  6f  is  arbitrary,  the  vector  in  the  parenthesis  must  be  zero,  or 
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[T  WT  + 


aMV 


[T]1||2 


-]1 


aQ 


l|[T]1||2  [T*Wj|o 
l|[T]?||2  -  aQ 


(52) 


Substituting  (38)  with  ^  defined  by  (35)  into  (52),  we  obtain 


l  a. [T  WT  +  0V]* 
1-1 


*v*  ■> 

y[t  w)g 


o 


(53) 


where 


IIIT]1|P  -  aQ 


(54) 


l|[T]f  |p 

Y - - - - - -  (55) 

l|[T)f |p  -  aQ 


Premultiplying  (53)  by  ^  and  using  the  orthogonality  relationships  (36)  and 
(37)  we  obtain 


a1(A1  +  6)  -  YC1 


(56) 


where  is  defined  by  (42).  Substituting  these  into  (38),  we  have 


I 


N 

l 

i-1 


Ai  + 


6  i 


(57) 


Next,  using  (57)  and  the  orthogonality  relationships  (36)  and  (37),  we  can 
write  the  Q  defined  by  (15)  in  terms  0  as 


imp 

IlmllP 


N 

l 

i-1 


(*i  + 


6) " 


N 

l 

i-1 


IK 


1  "i 

+  B)2 


(58) 


With  (57),  the  pattern  synthesis  error  e  extracted  from  (50)  is  given  by 
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2 


N 


N  |C±| 2 


»?  '“il 


'  n  Jj  <»,  +  6)-  "  Y  Jj  »,  +  6  "  v1i1»1  *  8  +  Hgi 


||2  (59) 


Manipulation  with  (54),  (55),  and  (58)  results  In 


?  All 

i-i  *i  +  8 


N 


(60) 


Equation  (60)  Is  more  easily  obtained  by  setting  to  zero  the  partial 
derivative  of  (59)  with  respect  to  y  .  Substituting  (60)  Into  (59), 
we  finally  have 


e 


'cii2  V 

(xt  +  S)  ) 

353” 

(X4  +  8) 2 


(61) 


It  is  desired  to  find  a  >  0  such  that  (48)  is  satisfied.  Since  (58) 
expresses  Q  not  in  terms  of  a  but  in  terms  of  8,  a  relation  between  a  and 
6  must  be  found.  Solving  (54)  for  a,  replacing  Q  by  expression  (58),  and 
recalling  that  || Tf  ||  is  proportional  to  Y  given  by  (60),  we  obtain 


a 


i  ?  iii!L 

M  if1  \±  +  t 


(62) 


The  derivative  of  (62)  given  by 


da 

d8 


I  f  _ 

Mi-l^i 


lcil: 


i 

w 


(63) 


indicates  that  the  right  hand  side  of  (62)  is  a  monotone  increasing  func¬ 
tion  of  8  except  at  6  -  -  X^  for  i-l,2,...N  where  it  jumps  suddenly  from 
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+“  to  -®.  If  the  eigenvalues  X^  are  arranged  in  decreasing  order,  then 
for  a  fixed  a  >  0,  (62)  has  one  root  6  outside  the  interval  (-X^,  0). 

There  is  also  a  root  in  each  of  the  intervals  (-X^,  -X1+^)  for  i-1,2, . . .N-l. 
Each  of  these  N  roots  represents  a  stationary  point  of  the  Lagranglan  J. 
Substituting  (58),  (61),  and  (62)  into  (47),  we  obtain 

J-lli0ll2-y  (64) 

This  shows  that  the  latter  N-l  roots  in  the  intervals  (-X^,  -X  )  are 
extraneous  because  tiey  all  render  J  larger  than  its  value  at  the  first 
root  outside  the  interval  (-X^,  0). 

We  seek  0  outside  the  interval  (-X^,0)  such  that  expression  (58)  is 
equal  to  Qq.  Since  Q  is  a  monotone  decreasing  function  of  a  and  a  is  a 
monotone  increasing  function  of  0  outside  the  interval  (-X^,  0),  it  follows 
that  Q  is  a  monotone  decreasing  function  of  0  outside  the  interval  (-X  ,  0) . 
Hence  there  is  only  one  0  outside  the  interval  (-X^,  0)  for  which  expression 
(58)  is  equal  to  Qq. 

That  Q  Is  a  monotone  decreasing  function  of  0  outside  the  interval 
(-X^,  0)  can  be  shown  from  (58)  as  follows: 


42 

dB 


-2M 


N  ICJ2  N  |C1|2X1  N  |C1|2X1  N  ICJ2 

1-1""(x1+B)Y  i-i  “  i-i  i-i  TV557 


N 

(  l 


i-l<V8) 


4) 2 


(65) 


Replacing  the  product  of  sums  in  the  numerator  of  (65)  by  a  double  sum, 
we  can  write  (65)  as 


42  .  _ 

dB 


2M 


N  N  |Ci|2|C1|2(Xi-X1) 

A  +  6)1 

- n  Tc/K. - - 

(I 


i-l<*i  +  8) 


?)  ‘ 


Combining  the  (i,j)  and  (j,i)  terms,  we  have 


(66) 
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(67) 


£2. 

dB 


-  2M 


?  J;1  |C1!2|C1I2(X1  - 

1-2  i-i  <y+  8?^  ♦  S)1 


N 

(  I 


lCi^i 


i-1  (Xi 


+  0)‘ 


Expression  (67)  is  negative  when  6  is  outside  the  interval  (-A^,  -A  ) , 
hence  Q  is  a  monotone  decreasing  function  of  B  outside  the  interval 

(-X1»  "Xn^* 

If  Qq  is  larger  than  the  Q  of  the  source  function  which  minimizes  c 
when  there  is  no  constraint,  namely  if 


Q 


o 


>  M 


N 

l 


1-1 

N 

l 

i-1 


(68) 


then  the  optimum  source  function  is  the  unconstrained  optimum  f  obtained 
by  setting  6  -  0  in  (57)  and  (60).  However,  if  Qq  is  less  than  expression 
(58)  at  a  *  »  corresponding  to  6  -  -  A^,  namely  if 


Q0  < 
o  \1 


(69) 


then  it  is  impossible  to  satisfy  the  constraint  (46).  If 

N  |Cj|! 

L  x 

i-1  i 


(70) 


then  the  optimum  ?  is  given  by  (57)  and  (60)  where  B  is  the  unique  number 
outside  the  interval  (-A^,  0)  for  which  expression  (58)  is  equal  to  Qq.  This 
B  is  easily  computed  using  Newton's  method  with  the  understanding  that  the 
iterations  are  not  allowed  to  proceed  into  the  forbidden  interval  (-A^,  0).  If 
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(71) 


r  1  Ql  M 

—  O  — 


N 

I  IcJ2 

1-1  1 

l 

1-1  1  1 


then  8  Is  negative,  but  if 
N 

l  IcJ2 

..  1-1 


l  lCJ\ 

1-1 


<  Q  <  M 
o 


N  |CtP 
i-lT” 

i-i  i 


(72) 


then  8  is  positive. 

The  same  examples  used  to  Illustrate  the  constrained  j|f||2  solution 
were  run  to  Illustrate  the  constrained  Q  solution.  To  check  on  the  accuracy 
of  the  computer  program,  the  constraint  on  Q  was  chosen  to  be  the  Q  ob¬ 
tained  for  each  case  In  Table  6.  Hence,  the  constraint  was  Q  -  4.05  for 
case  (a),  Q  -  4.09  for  case  (b),  Q  -  3.55  for  case  (c) ,  and  Q  -  16.15  for 
case  (d).  The  final  synthesized  patterns  were  indistinguishable  from  those 
of  Fig.  9.  The  final  values  of  ||f||2  and  E  are  given  in  Table  7.  Note  that 
the  error  E  of  Table  7  is  always  less  than  or  equal  to  that  of  Table  6. 

Table  7.  Source  norm  squared,  normalized  pattern  error,  and  quality 

factor  for  pattern  synthesis  with  constrained  quality  factor. 


llfll2 

E 

Q 

coordinate  origin  a 

4.23 

0.324 

4.05 

coordinate  origin  b 

4.26 

0.325 

4.09 

coordinate  origin  c 

4.19 

0.234 

3.55 

field  phase  alternating 

1.08 

0.957 

16.15 

magnitude  synthesis 

8.51 

0.190 

6.82 
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VII.  DISCUSSION 


Only  a  few  examples  are  given  In  this  report,  hence  it  is  difficult  to 
draw  general  conclusions  on  the  synthesis  procedures.  However,  for  the 
examples  chosen,  the  normalized  pattern  synthesis  error  was  greater  than  0.4 
when  the  sources  were  X/2  or  X  apart.  This  indicates  that  they  cannot  radiate 
a  pattern  very  close  to  the  chosen  cosecant  pattern.  Even  when  thfe  sources  were 
X/4  apart,  the  normalized  pattern  error  was  of  the  order  of  0.3  when  the  phase 
of  the  radiated  field  was  specified.  The  field  magnitude  pattern  synthesis 
procedure  reduced  this  error  to  C.172.  Constraining  the  source  norm  squared  or 
the  quality  factor  to  about  1/3  of  its  unconstrained  value  increased  this  error 
to  only  0.191.  Hence,  it  is  possible  to  reduce  ||f||2  or  Q  by  factors  greater 
than  3  with  little  change  in  the  synthesized  pattern. 

The  unconstrained  least-squares  pattern  synthesis  procedure  gives  the 
source  vector  (16).  When  this  is  substituted  into  (23),  we  obtain  for  the 
pattern  synthesis  error 


e  -II1J2  -|I(T]?I|2  (73) 

If  the  pattern  error  is  large,  the  synthesized  pattern  [T]f  must  have  a  small 
norm  compared  to  that  of  g^.  The  worst  possible  case  would  be  that  for  which 
the  space  of  [Tj?  is  orthogonal  to  gQ.  In  this  case  the  minimum  synthesis 
error  would  be  ||g0l|2  and  the  source  vector  would  be  1  »  0. 

The  constrained  norm  pattern  synthesis  gives 

-  N  l°J% 

II  miW  ■  (74) 

which  is  certainly  less  than  the  norm  squared  of  the  pattern  obtained  from 
unconstrained  pattern  synthesis.  If  the  constrained  Q  source  vector  is 
written  as  Y?,  then  Y  must  be  chosen  to  minimize 

e  -  II  Y  (T]f  -  gj2  (75) 
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Setting  to  zero  —y  we  obtnin 


W]gQ 

ymimr 

c(Y)  -  Hill2  -  ||y{t]?!|2  >  o 


(76) 

(77) 


so  that  the  norm  squared  of  the  constrained  Q  synthesized  pattern  is  less 
than  or  equal  to  ||g0ll2* 

The  field  magnitude  synthesis  procedure  of  section  IV  can  be  looked 

at  from  a  different  point  of  view.  Expression  (16)  implicitly  defines  an 

operator,  say  P,  which  gives  the  phases  of  the  elements  of  1  in  terms  of 

the  phases  0  of  the  elements  of  g  when  the  magnitudes  h  of  the  elements 
in  o  in 

of  gg  are  fixed.  The  iterative  procedure  of  section  IV  successively  oper¬ 
ates  with  P  on  an  initial  vector  of  phases. 


Convergence  is  obtained  when 


\  -  p<V 

t2  -  P(11) 


n+1 


n+1 


■  Pl\> 

approaches 


Evidently, 


lim  S 

n-H» 


n 


(78) 


(79) 


is  a  fixed  point  of  the  operator  P.  The  successive  operations  with  P 
could  have  been  replaced  by  Newton's  method  in  which  the  previous  iterate 
is  improved  by  adding  AS 
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(80) 


P(2  +  46)  -  6  +46 

n  n 

-+•  AP  “►  -► 

p(6  )  +  ~  46  -  6  +46 


(il_  _  i)A8  -  1  -  p(|  ) 

86  n  n 

n 

Presumably,  Newton's  method  requires  fewer  iterations,  but  each  iteration 

3P 

is  complicated  by  the  expression  ( -  1)  which  is  actually  a  square  matrix 

36 

to  be  computed  and  inverted.  Pn 

The  eigenvectors  $  of  (35)  will  be  real  if  the  elements  of  [¥*T]  are 
real.  The  following  development  reveals  the  circumstances  under  which  the 
elements  of  ft  T]  are  real,  or  nearly  real.  Using  (18), 


(T*  T)  -  I  e  mn  10 

mn  1.1 


where 


*  ^(*  ~  *)  +  (y  -  y„)‘ 

mn  m  n  m  n 


-1  yn  ~  ym 

,  i  n  m 

4>  *  tan  - 

o  x  -  X 

n  m 


Using  the  wave  transformation  [11] 


jkp _ eos 

mn  i  o 


r  a  JqOK-*  ) 

l  jq  J  (kp  )e 
q  mn 
q«-«  ^ 


one  obtains 


cc  -Jq<p„  M  jq4>. 

I  jq  J  (kp  )r  l  e 

q  mn  ,, 

M  i»l 


[11]  R.  F.  Harrington,  "Time-Harmonic  Electromagnetic  Fields,"  McGraw-Hill 
1961,  p.  231,  Eq.  (5-101). 
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Since  the  ere  equally  spaced  and  because  can  be  absorbed  Into  $qI 
we  are  at  liberty  to  take 


In  which  case 


♦l  '  (1-»  f 


M  jq4>.  Ch 

l  e  1  -■ 

1-1  0 
v_ 


q  -  an  Integer  multiple  of  M 
otherwise 


(86) 


(87) 


whence 

r  Mo  "jMq*0 

(T  T)  -Ml  j"q  JM  (kp  )e  0  (88) 

mn  nq  mn 

q— oo 

or 

00 

(T* T)  -  M[J  (kp  )  +  2  l  jMqJ  (kp  )cos(Mqi  )]  (89) 

mn  o  mn  ,  nq  mn  o 

q-l  M 


-v*  a/* 

Hence  (T  T)  is  real  whenever  M  is  even.  If  M  is  odd,  then  (T  T)  is 
mn  mn 

nearly  real  whenever  M  >>  kp  because  from  the  table  on  page  A07  of  [12] 

mn 

Jn(x)  is  very  small  when  n  >>  x. 

Newton's  method  was  used  to  find  tne  root  a  of  (45)  and  the  root  8  of 
(48)  where  Q  Is  given  by  (58).  However,  if  the  starting  value  Is  far  from 
the  root,  the  first  few  Iterates  of  Newton's  method  are  probably  not  very 
well  directed  because  the  change  In  the  variable  a  or  8  from  one  iteration 
to  the  next  Is  not  a  good  indicator.  In  the  beginning,  an  Interval  halving 
procedure  would  probably  give  as  much  improvement  per  iteration  and  take  a 
lot  less  time  per  iteration  because  the  derivative  need  not  be  calculated. 
However,  the  interval  halving  procedure  requires  a  negative  and  a  positive 
value  of  the  function  in  order  to  start.  Also,  it  would  be  difficult  to 
decide  when  to  change  from  interval  halving  to  Newton's  method.  Newton's 


[12]  M.  Abramowltz  and  I.  A.  Stegun,  "Handbook  of  Mathematical  Functions  with 
Formulas,  Graphs,  and  Mathematical  Tables,"  National  Bureau  of  Standards, 
1964,  pp.  355-433. 
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method  may  be  modified  by  replacing  the  derivative  by  a  finite  difference 
approximation,  but  then  it  is  feared  that  the  time  saved  per  iteration  may 
be  offset  by  slower  convergence. 

The  computer  programs,  with  operating  instructions  and  sample  input- 
output  data,  for  all  the  examples  of  this  report  will  be  given  in  Scientific 
Report  No.  3  of  this  contract.  It  is  hoped  that  further  examples  wl/11  be  run 
in  the  future  to  better  assess  the  capabilities  of  these  synthesis  programs. 
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APPENDIX 


ALTERNATIVE  METHOD  FOR  FIELD  MAGNITUDE  PATTERN  SYNTHESIS 


An  alternative  method  for  field  magnitude  synthesis  is  as  follows. 
Expression  (21)  can  be  rewritten  as 

e(f,B)  -  -  2  Real(e  w^  [  f  T  ) 

1  n-1 


N 

2  M 

l  f  T. 

S  n  in 
n-1 

+  I  w 
B=1  " 

m^i 


The  minimum  of  expression  (A-l)  with  respect  to  8^  occurs  when  is  the 
N 

angle  of  J  f  T^  .  Thus  *8  possible  to  minimize  (21)  with  respect 
n-1 

to  Bj  for  i  «  1,2,... M.  The  proposed  alternative  method  for  field  magnitude 
synthesis  consists  of  minimizing  (21)  successively  with  respect  to 
8^,  ^’''‘^M’  8^,  82»>«-8mi  etc.  Hence,  one  angle  is  changed  at  a  time  in 
this  method,  compared  to  all  angles  being  changed  at  once  in  the  method  used 
in  the  text.  It  is  not  known  which  method  converges  faster,  since  tests  were 
not  made. 
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